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Excitation transfer through interacting systems plays an important role in many areas of physics, chemistry, 
and biology. The uncontrollable interaction of the transmission network with a noisy environment is usually 
assumed to deteriorate its transport capacity, especially so when the system is fundamentally quantum mechani- 
cal. Here we identify key mechanisms through which noise such as dephasing, perhaps counter intuitively, may 
actually aid transport through a dissipative network by opening up additional pathways for excitation transfer. 
We show that these are processes that lead to the inhibition of destructive interference and exploitation of line 
broadening effects. We illustrate how these mechanisms operate on a fully connected network by developing a 
powerful analytical technique that identifies the invariant (excitation trapping) subspaces of a given Hamiltonian. 
Finally, we show how these principles can explain the remarkable efficiency and robustness of excitation energy 
transfer from the light-harvesting chlorosomes to the bacterial reaction center in photosynthetic complexes and 
present a numerical analysis of excitation transport across the Fenna-Matthew-Olson (FMO) complex together 
with a brief analysis of its entanglement properties. Our results show that, in general, it is the careful interplay 
of quantum mechanical features and the unavoidable environmental noise that will lead to an optimal system 
performance. 



Transport phenomena have been central to quantum me- 
chanics since its early days [1]. This pivotal role has been 
recently renewed by the prospect of transferring quantum in- 
formation across quantum networks [2] and the recurring in- 
terest in understanding the fundamental processes influencing 
energy transport in photosynthetic molecules Light har- 
vesting complexes can harness the available light energy at ef- 
ficiencies well above 90% in conditions that are often defined 
as hot and wet [4, 5]. Understanding the workings of such 
a process has the potential to be immensely valuable from a 
technological point of view, potentially paving the way for the 
design of novel nanofabricated structures for quantum trans- 
port and optimized solar cells. There is already an impor- 
tant body of work analyzing energy transport in these systems 
[6, 7, 8, 9, 10, 11]. However, the clear identification of the 
basic mechanisms underlying this remarkably efficient natural 
transfer scheme has remained elusive. Recently, a sequence of 
beautiful experiments using very effective nonlinear spectro- 
scopic techniques probed the dynamics of delocalized exciton 
states in light-harvesting complexes [12], the Fenna-Matthew- 
Olson (FMO) complex [13, 14] and conjugated polymer sam- 
ples [15], the latter modeling multichromophoric systems. 
Based on these experimental observations, quantum coher- 
ence across multiple chromophoric sites has been suggested 
as the probable cause of the highly efficient energy transfer in 
photosynthetic systems. The influence of initial exciton de- 
localization in the efficiency of excitation transfer was inves- 
tigated by [ 1 6] and it has been suggested that the enhanced 
transfer rates may be attributed to the exploitation of princi- 
ples of quantum search algorithms by the quantum dynam- 



ics of the FMO complex [3, 14]. At the same time, how- 
ever, excitations in chromophoric complexes are subject to 
strong dephasing noise from a quasi-continuum of environ- 
mental phonon modes, which, in turn, suggests that quantum 
coherence will be short-ranged rather than extending across 
the entire complex. This short ranged nature of quantum cor- 
relations may suggest that concepts from quantum computa- 
tion will not play as decisive a role as it is generally accepted. 
In fact, efficient pure state quantum computation requires long 
range entanglement [17]. It therefore remains unclear whether 
the presence of quantum entanglement has any specific func- 
tional role in the process of energy transport across biolog- 
ical systems, or it can simply be viewed as the unavoidable 
by-product of the existence of a coherent quantum evolution 
during a certain time scale in such systems. 

Inspired by these experimental results and their interpreta- 
tions, further theoretical studies of energy transport in light 
harvesting complexes have been carried out recently [18, 19]. 
They investigate the role of noise, and in particular dephas- 
ing, in the process of exciton transport in these complexes. 
Contrary to the conventional wisdom that noise always de- 
teriorates the performance of a system, several instances are 
known, in classical as well as quantum mechanics, where 
noise is known to be of advantage in enhancing the system's 
response. In the classical phenomenon of stochastic resonance 
[20], a non-linear system may benefit from the presence of 
noise to achieve an enhanced sensitivity to a weak signal. Mo- 
tivated by this observation, it was successfully demonstrated 
that noise may have the capacity to generate genuinely quan- 
tum mechanical features such as quantum coherence and en- 
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tanglement [21, 22, 23]. Taking this further, it was shown re- 
cently that noise, in the form of local dephasing, has the abil- 
ity to enhance the rate and efficiency of energy transfer when 
compared to a perfectly quantum coherent system [18, 19]. 
Indeed, the highly effective exciton transfer in light harvest- 
ing complexes requires the simultaneous presence of quantum 
coherent evolution and dephasing noise. Simulations have 
also shown that, in the competition between coherent dynam- 
ics and incoherent energy dissipation (population decay), the 
former regime lasts for about several hundreds of fs [24] at 
the beginning of the exciton transfer process that spans about 
5 ps [ I I ]. As a result, in addition to the expectation that the 
quantum coherence is short-ranged, it is in fact short-lived as 
well, and so will be any entanglement across the systems. 

In spite of these qualitative and quantitative successes based 
on numerical work, the underlying mechanisms by which 
noise supports transport have yet to be precisely identified and 
exemplified in simple models. Achieving this may assist us in 
replicating in artificial structures what nature does so well - 
energy transfer in light-harvesting complexes. 

In this paper we identify two fundamental mechanisms un- 
derlying dephasing-assisted transport and elucidate them in 
the context of two basic models. We will show that these are 
firstly, the suppression of destructive interference by adding 
local decohering noise such as dephasing or static disorder, 
and secondly, the enhancement of transfer by line broadening 
caused, for example, by the fluctuation of energy levels. Both 
mechanisms have the principal effect of opening additional 
channels for transport in the system; channels that would be 
inhibited under solely coherent evolution. 

We start by presenting in section I an abstract network 
model for quantum transport where local sites are subject to 
both local dissipation and dephasing and the excitation trans- 
fer to a reaction center is modeled via an irreversible cou- 
pling to a privileged trapping site, an excitation sink. In sec- 
tion II we introduce a simplified model in terms of a fully 
connected network, where all sites are coupled to each other 
with equal coupling strength. This model will allow us to 
develop a powerful analytical technique to elucidate whether 
a given Hamiltonian is susceptible to support noise assisted 
transport. This will be done by characterizing the invariant 
(excitation trapping) subspaces of the Hamiltonian. With the 
aim of making the pace of the paper reasonably fluid, the de- 
tailed procedure is presented in a separate Appendix where 
we develop a technique for solving the complete master equa- 
tion for the fully connected network. Using these analytical 
results, we show in sections II.A-II.D that dephasing assists 
transport in fully connected networks while dissipation only 
does so in some cases. We provide examples when pure dis- 
sipation may assist the transport through a network. Section 
lI.E revises line broadening as a classical effect. In section III 
we extend the model to account for forms of spatially corre- 
lated noise. This theoretical analysis allows for the identifi- 



cation of the fundamental mechanisms underlying dephasing- 
assisted transport which are then applied in section IV to the 
specific example of energy transfer across a light harvesting 
system, the FMO complex, which may be modelled as a net- 
work of seven nodes, fully connected, albeit with non-uniform 
coupling strengths. We show how experimental results con- 
cerning exciton transfer cannot be reproduced in the absence 
of dephasing while the inclusion of this form of local noise 
boosts energy transfer towards the observed values in the cor- 
rect time scale. We show that the two mechanisms presented 
here by which dephasing assists excitation transfer are very 
robust against variations of the detailed structure of the sys- 
tem. The theoretically expected and numerically observed 
weak dependence on the noise strength suggests that exper- 
imental results which have been obtained at 77 K should re- 
main broadly valid at higher temperature [14] and that these 
processes will be robust against perturbations due to different 
environments. We conclude section IV with an evaluation of 
entanglement generation and transmission within this simple 
noise model, including local and correlated dephasing, and 
show that quantum correlations are expected to be confined to 
the initial steps of the transport process and die out well before 
the excitation is fully transferred. In section V, we summarize 
the main physical ideas behind our approach in an intuitive 
form and we conclude in section VI by suggesting that the 
principles outlined here may find broad applicability in under- 
standing the dynamics of biological complexes and transport 
phenomena and discuss future work in this direction. 

I. THE NETWORK MODEL 

Light-harvesting complexes are typically constituted of 
multiple chromophores which transform photons into exci- 
tons and transport them to a reaction center [25]. Experi- 
mental studies of the exciton dynamics in such systems re- 
veal rich transport dynamics consisting of short-time coherent 
quantum dynamics which evolve, in the presence of noise into 
an incoherent population transport which irreversibly trans- 
fers excitations to the reaction center In order to elucidate the 
basic phenomena clearly without overburdening the descrip- 
tion with detail, we consider the relevant complexes as sys- 
tems composed of several distinct sites, one of which is con- 
nected to the chromosomes while another is connected to the 
reaction center This complex effective dynamics will then be 
modelled by a combination of simple Hamiltonian dynamics 
which describe the coherent exchange of excitations between 
sites, and local Lindblad terms that take into account the de- 
phasing and dissipation caused by the external environment. 
A network of N sites will be described by the Hamiltonian, 

N 

H = Y^ /it^jCr+cr- + hVj^i{(J-(J+ + cr+crf), (1) 
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where ('^J^ ^^e raising and lowering operators for 
site j, hctij is the local site excitation energy and v^.i de- 
notes the hopping rate of an excitation between the sites k 
and /. We will also designate a site representing the zero 
exciton state of the complex, which appears in operators such 
as aj' = \j) (0|, where the state \j) denotes one excitation 
in the site j. We will assume that the system is susceptible 
simultaneously to two distinct types of noise processes, a dis- 
sipative process that transfers the excitation energy in site j to 
the environment (with rate Tj) and a pure dephasing process 
(with rate 7j) that destroys the phase coherence of any super- 
position state in the system, i.e randomizing the local excita- 
tion phase in site j. Following previous studies on dephasing- 
assisted transport [18, 19], we describe both processes using 
a Markovian master equation with local dephasing and dissi- 
pation terms [26]. As we shall see, the mechanisms that we 
discuss here are quite insensitive to the details of the noise 
model. In the Markovian master equation approach, the dis- 
sipative and the energy-conserving dephasing processes are 
captured, respectively, by the Lindblad super-operators 

N 

C-,uss{p) = ^rj[-{a+a-,p} + 2a-p4], (2) 

N 

Finally, the total transfer of excitation is measured by the pop- 
ulation in the 'sink', numbered + 1, which is populated by 
an irreversible decay process (with rate Fjv+i) from a chosen 
site k as described by the Lindblad operator 

{o"i^^jv+iO"w+i'^fe .P}]- (4) 

For definitiveness and simplicity, the initial state of the net- 
work at t = will be assumed to be a single excitation in 
site 1 (i.e., state |1)), unless stated otherwise. The model 
is completed by introducing the quantity by which we mea- 
sure the efficiency of network's transport properties, that is 
the population transferred to the sink Psink{t), which is given 
bypsmfe(i) = 2rjv+i Jq pkk{t')dt' . 

II. THE FULLY CONNECTED NETWORK (FCN) 

The FMO complex can be described by a Hamiltonian in 
which every site in the network is coupled to every other sites 
but with site dependent coupling strength [11, 19]. The domi- 
nant couplings are the nearest neighbour terms, but significant 
hopping matrix elements exist also between more distant sites. 
This motivates the study of networks with a high level of con- 
nectivity (see also Ref. [27] for quantum state transfer), and 
in this section, for simplicity, we will look in detail at a fully 
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FIG. 1: Paink VS. time is shown for a fully connected graph of 
N = 10 nodes with F^ = for i = 1, . . . Af, J = 1, Fjv+i = 1. At 
t = one excitation is in the site 1 and the energy transfer evolves ac- 
cording to the quantum dynamics investigated in the text. The cases 
of an energy mismatch, i.e. uji = 1, cji^ti = (red line), a dephas- 
ing mismatch, i.e. 71 = 1, ^i^i = (green line), and the basic case 
= iOi = (black line) Vi are shown. For the latter, the popula- 
tion in the sink asymptotically reaches psink = 1/9. Note that both 
the red and the green lines reach unit transfer asymptotically, but at 
considerably different rates. 



connected network (FCN). The FCN is characterised by equal 
hopping strengths between all sites, i.e. hvjj = J for any 
j 7^ /. Remarkably, for the case of a uniform network, i.e. 
one in which uij, -fj, and Tj are the same on every site, an ex- 
act analytical solution can be found for the density matrix of 
arbitrarily large networks. These exact solutions are obtained 
by defining a set of collective variables, and the formal devel- 
opment of the solutions is given in the Appendix. These exact 
solutions provide insight into the real-time transport dynamics 
of the FCN, and, by considering different scenarios in which 
various noise effects are present or absent, we can isolate the 
various mechanisms that contribute to the dephasing-assisted 
transport, as discussed below. Moreover, having isolated these 
mechanisms, we find that they can be described in an intu- 
itive wave mechanics picture which highlights the relevance 
of these processes in a much wider range of quantum network 
models. 



A. No dissipation, no dephasing: Destructive interference 

This case is described by -fj = F^ = and, for simplicity, 
ujj = for any j ~ 1, . . . ,N. The only irreversible process 
left in this system is the decay of population to the sink from 
site A^. The exact solution (in the Appendix) predicts the strik- 
ing result that as t ^ 00 the total amount of the excitation that 
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is transferred to the sink is given by 



Psinkioo) 



1 



N - 1 



(5) 



In Fig. (1), the time evolution of the sink population is shown 
(black hne), for the case of = 10. For any reasonably large 
network the transfer to the sink is very small, even in the limit 
t ^ oo. This should be contrasted with classical hopping, 
e.g. a random walk, model, in which the excitation can be 
shown to be completely transferred to the sink as i oo. 
The difference between these results is a consequence of the 
wave-like nature of the quantum dynamics in the FCN and 
the presence of destructive interference. This can be seen by 
considering the final (i.e., t — > oo) state of the network given 
by 



Pit 



where 



N -1 



1 



N - 1 



|iV+l)(7V + l| (6) 



I*) 



V{N-l){N-2) 



(7) 



and |1) is the initially occupied site. Even though each indi- 
vidual site has a finite amplitude J for transfer to site A^, the 
amplitudes coming from |1) cancel those from |j) due to the 
minus sign in the superposition. As a result of this destruc- 
tive interference, the net amplitude connecting j^E') and \N) 
must be zero, which means that the excitation stored in j^E') 
is unaffected by any process that acts locally at N. Thus a 
quantum FCN can protectively store some of the excitation 
in superpositions that have no effective overlap with site N, 
as illustrated in Fig. (2(a)). A powerful way of understand- 



(a) 



(b) 



Sink 

FIG. 2: A fully-connected three-site network. In (a) the excitation 
wave function is delocalised over two sites (red and green) with 
equal probability of being found at either site. However, as the wave- 
function is antisymmetric with respect to the interchange of red and 
green, this state has no overlap with the dissipative site (blue) due to 
the destructive interference of the tunnelling amplitudes from each 
site in the superposition. The network can therefore store an excita- 
tion in this state indefinitely. In (b) pure dephasing causes the loss 
of this phase coherence and the two amplitudes no longer cancel, 
leading to total transfer of the excitation to the sink. 



ing how this final state emerges without actually solving the 
full dynamics is to invoke the notion of an invariant subspace. 
We define an invariant subspace to consist of a set of states 
which are both eigenstates of H and have no overlap with site 
N. The eigenstates of the FCN consist of a state of the form 
^N-i J2f=i \j) with energy E = J{N - 1) and iV - 1 
degenerate eigenstates with E = —J. Due to the degeneracy 
of these states, they can be represented in several ways. One 
particularly useful representation consists of anti-symmetric 
superpositions of pairs of sites of the form \ipj) = |1) — |j) 
where j = 2 . . . N. We now expand the initial system state 
in terms of the invariant subspace and the remaining states to 
find 



ID = 



N - 1 



Invariant 



' N-1 

E(l^.> I +VN\cj^)-\N) 

J=2 



(8) 



The subsequent time evolution of the system is very sim- 
ple. The states in the invariant subspace are eigenstates of 
the Hamiltonian and are not affected by the open-system dy- 
namics that only acts at site A'^. Therefore, their evolution is 
purely coherent and, being degenerate, can be described by a 
simple global phase. The remaining components of the ini- 
tial state expansion evolve in the non-invariant subspace de- 
fined by having a finite overlap with site N and are therefore 
transferred into the sink. The net result is that at long times 
the weight held in the network is simply that contained in the 
invariant part of the initial state expansion. A further conse- 
quence of this analysis is the prediction that, if the system is 
initially prepared in one of the invariant states, then there are 
no dynamics: the excitation is completely trapped in a noise- 
protected stationary state and Psink{oo) = 0, as illustrated in 
Fig. (2(a)). Note that this phenomenon is reminiscent of co- 
herent population trapping in quantum optical [2S] and con- 
densed matter [29, 30] systems. 

The power of this approach is that the argument given above 
must also be valid for other Hamiltonians provided they con- 
tain an invariant subspace. Invariant subspaces turn out to be 
rather easy to find, and can be systematically generated for 
any network with some degeneracy. This would naturally im- 
ply that large invariant subspaces can be found in systems of 
high symmetry, and the larger the invariant subspace is, the 
larger is the spectral weight retained by the network. A brief 
discussion of more general networks within the invariant sub- 
space approach can be found in the Appendix. These results 
are independent of the noise model provided that it only acts 
locally at N . However, under non-local noise it may also be 
possible to find invariant subspaces as well. 
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B. Static disorder: Suppressing destructive interference 

Another advantage of the invariant subspace approach is 
that it can be extended to situations in which the uniformity 
of the network is perturbed and allows us to study one of 
the key mechanisms for suppressing destructive interference. 
Changes to the local site energies are local perturbations, and 
it should once again be possible to construct an invariant sub- 
space that does not feel their influence. This requires that we 
define a new invariant subspace that is spanned by all eigen- 
states that do not have any overlap with either site N, con- 
nected to the sink, or the perturbed sites. For the FCN this 
amounts to a reduction of the previous invariant subspace by 
the number of perturbed sites, that is, we now exclude all \ 
corresponding to the perturbed sites j. The final state of the 
system now consists of the part of the initial state expansion 
that lies in the reduced invariant subspace, which for FCN im- 
mediately implies 

where D is the number of sites with different site energies 
[31]. This implies that a network where all the site energies 
differ {D ^ N — 2) results in Psink{oo) ~ 1, i.e. complete 
excitation transfer. A natural consequence of this result is that 
a network with randomly disordered energies always leads to 
perfect transfer as t — > oo. 

While in the FMO complex and biological complexes in 
general other mechanisms may also play a role, the above 
strongly suggests that it is not an accident that in many of 
those complexes a moderate amount of disorder is present in 
the distribution of the site energies. 

The significance of this observation should not be underes- 
timated. It contrasts with the idea that disorder leads to ei- 
ther weak or Anderson localization in the system [32, 33], 
and that static disorder always inhibits transfer. That this is 
not the case in the FCN is already demonstrated by the above 
results where at least a certain amount of disorder is essential 
for achieving complete population transfer The origin of this 
perhaps surprising behavior is that, unlike systems described 
by the Anderson tight-binding Hamiltonian, the FCN trans- 
port to the sink is already strongly suppressed in the absence 
of disorder due to destructive interference of transition am- 
plitudes to the state coupled to the sink. As we have shown, 
this destructive interference is inhibited by any finite amount 
of disorder and thus disorder can initially lead to an enhance- 
ment of transport in the FCN. However, while finite disorder 
leads to PsmA:(oo) = 1, the rate at which the excitation is 
transferred to the sink is not a monotonic function of the dis- 
order strength - see Fig. (3). For weak disorder the spectrum 
of the FCN still contains states that are fairly close to the un- 
coupled subspace, i.e. they couple only weakly to the sink, 
and transport through the network will be slow. As the size 



of the disorder increases, the eigenstates become increasingly 
distant from the uncoupled subspace and should lead to an in- 
crease in the transport rate. Increasing the magnitude of the 
disorder further eventually causes the transfer through the net- 
work to become slower again due to the reduction in overlap 
between the sites, as will be discussed later on. 



C. Local dephasing without dissipation: Suppressing 
destructive interference 

The previous sections showed that the noise-free FCN can 
trap a large part of the initial excitation within a manifold of 
states that do not overlap with site and hence do not lead to 
transfer to the sink. These states are coherent superpositions 
of excitons in different sites, and are invariant in time due to 
destructive interference of the transition amplitudes. Static 
disorder however may suppress this destructive interference, 
thus releasing all or part of the trapped excitons. Here we 
show that local dephasing noise has a very similar effect be- 
cause it affects the relative phase in the states that are invari- 
ant under the Hamiltonian dynamics, e.g. mapping the invari- 
ant state (|1) - |2))/V2 to the evolving state (|1) + |2))/V2 
(in the master equation picture presented above for example 
through the action of the operator (t^ctj^), and thus perturb- 
ing the important destructive interference. This observation is 
quite general and may be made for arbitrary dephasing mecha- 
nisms both due to classical or quantum environments. Indeed, 
in the presence of local dephasing on all sites the size of the 
invariant subspace discussed there vanishes and one finds that 
as t ^ oo, Psink tends to unity for any finite 7. This may be 
confirmed analytically for the case of a uniform local dephas- 
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FIG. 3: Dependence of Psink (t) at a fixed time t = 20 as a function 
of 7/rjv-i-i, in the case of = 5, J = 1, and lu = 0. The initially 
sharp rise is due to the increasing rapidity at which the invariant sub- 
space is destroyed, whilst the decreasing part is due to quantum Zeno 
effects. 
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FIG. 4: Population of the sink Psink (t) as a function of time t for in- 
creasing dephasing rates 7, in the case of Fiv-i-i = 1. For short times, 
the system dynamics are identical to the case of zero dephasing, and 
for very weak dephasing a significant fraction of the excitation can 
be trapped in invariant subspace for a time of approximately 7^^. 



ing rate — 7, for which the dynamical equations can again 
be solved exactly, as shown in the Appendix. This demon- 
strates the same fundamentally quantum mechanical aspect of 
dephasing-assisted transport as the Mach-Zender interferome- 
ter in Sec. V. A: the removal of destructive interference and the 
creation of new pathways for transport (see Fig. (2(b))). Given 
that all physical realisations of network models are subject to 
some dephasing, one might wonder why we have laboured 
to develop the notion of invariant subspaces, spaces which are 
anything but invariant in a realistic network. One reason is that 
the dynamics of the excitation transfer can still be strongly in- 
fluenced by the presence of destructive interference. This is 
shown in Fig. (4) which shows the evolution of Psink{t) for 
N = 5 and J = 1 as a function of time for various pure de- 
phasing rates. For 7 ^ Fjv+i the early time evolution is iden- 
tical to the case of no dephasing. This implies that the system 
evolves into the final state found in the last section and then 
slowly dephases leading to a gradual transfer of the remaining 
spectral weight to the bath. Interestingly, the speed at which 
the excitation is transferred is not a monotonic function of the 
dephasing rate. This is illustrated in Fig. (3) which shows 
Psink{t) as a function of 7/r7v+i at a fixed time of t = 20. 
The initial sharp increase in the transport rate is due to the de- 
struction of the invariant subspace as discussed above, whilst 
the slowing of the transfer rate for large dephasing rates can 
be ascribed to either a quantum Zeno effect or, equivalently, 
an effect of line broadening. We will discuss the latter in more 
detail below. 

In Fig. (1), we compare our base case (7 = cj = F = 0) to 
two other cases both of which lead to Psmk = 1- One is the 
case of just a local dephasing on the site of injection (7 = 
elsewhere), while in the other case we have a different site 



energy at the site of injection (w = elsewhere). We can 
see very clearly that pure dephasing leads to a faster transfer 
than pure energy mismatch. An insight of this nature could 
help us control the transport behaviour in laboratory and real 
systems much more effectively. As shown above by using the 
invariant subspace argument, one may increase the transfer to 
the sink in small controlled steps, thus carrying us to a limit 
where we can almost direct the flow of energy in systems at 
the microscopic level by engineering macroscopic parameters. 

D. Only local dissipation, no dephasing 

Dissipation of the exciton by the environment also leads to 
the suppression of destructive interference, which raises the 
natural question whether relaxation alone can enhance trans- 
port in the absence of pure dephasing or static disorder of site 
energies. The answer depends on the system and initial state. 
To address this question, we consider the case of jj ~ and 
Tj = F. The exact solution is presented in the Appendix and 
predicts 

J^Fjv+i(2r + Fat+i) 

- mTj,^,,j,N) ' ^^^^ 

where cj) is given by 0(F, Fat+i, J, N) = iT^ + SF^Fw+i + 

5T^Tl^, + rr%^, + r%^,J'{N - i) + t^j^n'- + 

TTn+i.PN^. From Eq. (10) we find ^P"g^^°°^ < and 
hence that pure relaxation always yields psinfe(oo) < {N — 
Here, the losses to the environment always offset any 
increase in the transport efficiency due to dephasing of the 
invariant states. For the case of both pure dephasing and re- 
laxation, the transfer is always sub-optimal as losses to the 
environment always lead to psi„fc(cx)) < 1. 

Interestingly, it should be noted that the conclusion that 
pure dissipation does not enhance transport to the sink is not 
always valid. To see this, consider a system as in Fig. (2(a)) 
that is initially prepared in one of the invariant states such as 
(|1) - |2))/\/2 so that one finds ps,„fc(oo) = for F = 0. 
Now let us assume that only site 2 suffers dissipation, F2 > 0. 
The time evolution is now composed of two contributions. In 
one part, the excitation is lost to the environment leading to 
no successful transport. The other part in which the exci- 
ton is not being lost is more interesting. This branch of the 
time evolution is described by a conditional Hamiltonian that 
is the original Hamiltonian supplemented, crucially, with the 
non-Hermitean additional term —iF2|2) (2| [34]. Then the ini- 
tial state ( 1 1) — 1 2) )/ V2 is not invariant under the conditional 
time evolution anymore and the relative weight of state 1 2) in- 
creases compared to state This in turn affects the destruc- 
tive interference for transitions to level |3) as the interfering 
amplitudes then do not have the same weight anymore. This 
releases the trapped population and leads to some transfer to 
site 1 3). 
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FIG. 5: Transport in a classical ladder like system of = 7 sites, 
with steps of unit size (Ai5i,i+i = 1), with a dissipation term (with 
rate Fjv+i) into a sink, connected to the site 7. The energy levels 
fluctuate with a Lorentzian distribution (with width 1) to take into 
account line broadening effects and a symmetric hopping rate given 
by + 1) = V(47r) l/(ASj,j+i2 + 1/16) between the sites j 

and i + 1 and Fat+i = 1 (red line), while also the noiseless case is 
shown (black line). Inset: Pdf distribution ioxpaink for t — 100 over 
a sample of 10* instances in both hopping scenarios. Notice that it is 
very narrow around the central Psink value. In both cases, it is clearly 
shown that line broadening dramatically enhances the transport. 



Thus in the absence of pure dephasing, it is possible to find 
instances of both positive and negative contributions to the 
transport as a result of relaxation. 

E. Line broadening as a classical effect 

So far we have placed an emphasis on the effect of suppres- 
sion of destructive interference, hence addressing the wave- 
like feature of the exciton transport. Dephasing may also be 
seen to broaden the resonance lines of the individual sites. 
This effect has a purely classical component as well. Indeed, 
fluctuating excitation energies of sites leads to stronger cou- 
pling between sites whose energy difference becomes small 
or even changes their energy ordering thanks to fluctuations 
in time [ I '>]. In order to clearly separate this mechanism from 
possible quantum interference effects discussed before, we ex- 
emplify it here using the purely classical model of N sites in 
an ascending ladder of sites as displayed in Fig (5) with the 
highest energy site connected to a sink. The probability of 
finding the excitation on the site i, i = 1, . . . , iV at a time t, 
Pi{t), is given by the Pauli master equation as 



dpijt) 
dt 

dPsinkjt) 

dt 



h,i+i Pi+iit) + h,i-iPi-i{t) (11) 

- {ki+i,t + ki-i,i + '^n+i5in) Pt{t) 
TN+lPNit) , (12) 



where fc^.j are the hopping rates from sites j into site i, re- 
spectively, r^r_|_i is the irreversible decay from the site N 
into the sink and we assume that the energy differences be- 
tween neighboring sites is AEi i^i = 1. We take account 
of line broadening by adding to this level spacing a noise 



term leading to a Lorentzian distribution of energy levels with 
width 1. Now, we consider two very closely related scenar- 
ios: namely (i) we consider thermal hopping, i.e. kij/kj^i = 
Q-^(Ei-Ej)/T ^jjjj j-jj-j jjjg hopping is symmetric ki j = kj i, 

with an hopping rate obeying a Lorentzian in the energy differ- 
ence of two neighboring sites. Notice that the latter is perhaps 
closer to the quantum case where the coupling between neigh- 
boring levels with an energy mismatch is better approximated 
by case (ii). In Fig. (5) the dramatic rise in transport due to line 
broadening is shown for the case of Lorentzian hopping rates. 
A similar behaviour was found also for other cases. Hence, 
we observe that energy mismatch in a system with static dis- 
order inhibits the transfer of an excitation up the ladder and 
the transport may be supported by addition of dephasing. 



III. CORRELATED NOISE 

The considerations so far have assumed Markovian local 
noise. In photosynthetic complexes that will be discussed in 
the following section it is however known that such a model 
is not fully accurate both due to the small size of the complex 
compared to the correlation length of the relevant bath and 
due to the fact that the time scales for interactions may be of 
the order of the bath correlation time. To gain an understand- 
ing of how spatial correlations affect decoherence, we study a 
Markovian master equation 



dp 
'dt 



(13) 



where the Hamiltonian and the dissipative term are chosen as 
Eq. (1) and Eq. (2) respectively while the Lindblad operator 
describing the correlated dephasing is given by 

Cdephip) = ^Jmn{~Crm<^m'^n<^nP- P'^n'^n'^m'^m 



'<^tn(^,nP<^n<^n + <^n <^n P(^m<^ra) ■ (14) 



Note that the Hermitian matrix of the 7„j„ must be positive 
semidefinite to describe a completely positive dynamics. The 
effect of the noise correlations can now be seen quite clearly 
by considering the special case of a two site system where 
711 = 722 = 7' 712 = 721 = a7, Lfe = for /s = 1, . . . , TV 
andH = 0. With the initial state p = i(|l) - |2))((1| - (2|), 
which led to destructive interference in Fig. (2), we find 

^ = 2(1 - a)7(|2) (2| + |2) (1| + |1) (2| + |1) (1|) . (15) 

For a = 1, i.e. perfectly correlated noise, we find p = 0, i.e. 
the state p remains stable. For anti-correlated noise a = — 1 
however this state would be decaying more rapidly than under 
uncorrelated noise. Hence, spatial correlations in the dephas- 
ing noise may act both to stabilize or destabilize specific co- 
herent superpositions and the related destructive interference 
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effects explained earlier in this work. Generally, of course 
the spatial correlations in the noise will not be perfect and de- 
pend for example on the distance between sites - more distant 
sites will tend to suffer weaker noise-correlations. Needless to 
say, the above master equation picture does not encompass all 
possible environments and also neglects non-Markovian ef- 
fects which may however be easily included for example by 
using time convolved master equations with a memory kernel 
that nevertheless preserve complete positivity [i^]. A thor- 
ough discussion of such effects goes well beyond the scope of 
this work as it would require us to take into account the details 
of the environment in our simulations. However, the present 
paper is aimed at clarifying fundamental mechanisms and we 
will therefore present a detailed study of the environment in a 
future publication. 

IV. LIGHT HARVESTING - THE FMO COMPLEX 
A. Excitation transfer 

Having elucidated the fundamental mechanisms that lead to 
enhanced transport in abstract, fully connected, networks, we 
will now consider them in the context of transport properties 
in the FMO complex [36] that has been the subject of several 
theoretical [18, 19, 24] and experimental [ ; ] studies recently. 
We will consider both local and non-local noise and will also 
compare the two settings to assess the impact of noise corre- 
lations. The FMO complex is a pigment-protein complex that 
funnels the excitation energy from the light-harvesting chloro- 
somes to the bacterial reaction center in green sulfur bacteria. 
It is a trimer of three identical units, each composed of seven 
chlorophyll a molecules embedded in a scaffolding of protein 
molecules. We model the FMO complex as a completely con- 
nected network, albeit with site dependent coupling strengths 
and site energies [11]. Specifically, the dynamics describing 
the FMO complex is composed of a Hamiltonian part describ- 
ing the coherent dynamics where the site energies and cou- 
pling constants have been taken from tables 2 and 4 of refer- 
ence [11]. We then find, in matrix form 



/ 


215 


-104.1 


5.1 


-4.3 


4.7 


-15.1 


-7.8\ 




-104.1 


220.0 


32.6 


7.1 


5.4 


8.3 


0.8 




5.1 


32.6 


0.0 


-46.8 


1.0 


-8.1 


5.1 




-4.3 


7.1 


-46.8 


125.0 


-70.7 


-14.7 


-61.5 




4.7 


5.4 


1.0 


-70.7 


450.0 


89.7 


-2.5 




-15.1 


8.3 


-8.1 


-14.7 


89.7 


330.0 


32.7 


V 


-7.8 


0.8 


5.1 


-61.5 


-2.5 


32.7 


280.0/ 



where the zero of energy has been shifted by 12230 for all 
sites, corresponding to a wavelength of = 800 nm (all num- 
bers are given in units of cm^^ =1.988865 • lO^^'' Nm = 
1.2414 10^** cV). The non-unitary part of the evolution is then 
described by Eq.(13). Recent work [ i , ] suggests that it is this 




time (ps) 

FIG. 6: Psink vs. time (in ps) for the FMO complex. We show 
the noiseless case (black line) and the transfer for optimal dephas- 
ing (red line). The optimal dephasing values are given by by 
{0.157, 9.432, 7.797, 9.432, 7.797, 0.922, 9.433} ps'^. The inset 
shows the very narrow probability distribution of transfer probabili- 
ties at t ~ 10 ps from 10^ samples, in the case of optimal dephasing 
but with a 20% static disorder in energy and coupling rates. This sug- 
gests that dephasing assisted enhanced transport is also robust against 
variations in the system parameters. 



site 3 that couples to the reaction centre at site 8. For this rate, 
somewhat arbitrarily, we chose g = 62.8/1.88 cm^^, cor- 
responding to about 6.283 ps^^ (note the equivalence h ~ 5.3 
cm~^ ps). The measured lifetime of excitons is of the or- 
der of 1 ns, which determines a dissipative decay rate of 
2Tk = 1/188 cm~^ and that we assume to be the same for 
each site [11]. 

As shown earlier, a completely coherent dynamics is often 
not most ideal for the transfer of excitons from the chloro- 
somes to the reaction center and does not match with exper- 
imentally observed transfer efficiencies. Indeed, the transfer 
in the FMO complex Pwsthecochloris aestuarii undergoing 
completely coherent dynamics caps at about 57%, shown in 
black in Fig. (6). The same figure shows the effect of local 
dephasing on all the sites. The enhanced transport is obtained 
with dephasing rates that were optimized numerically for a 
transfer time of about 5 ps [ i, i ']. The dynamical simula- 
tion of energy transfer in our model of the FMO complex is 
shown online through two movies in [37]. 

In the following we use the basic principles that we have 
formulated in the previous sections to understand where and 
how dephasing will be beneficial for the exciton transport in 
the FMO complex. These expectations will then be corrobo- 
rated by determining numerically the optimized local dephas- 
ing rates that will indeed be broadly in-line with the choices 
suggested below. When studying the dephasing-free evolu- 
tion, one observes a behaviour quite similar to that shown in 
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the black curve of Fig. (6). After a rapid rise in transfer prob- 
ability the system then enters a phase of slowed transfer that 
is due to the system having evolved into a state that is almost 
invariant in time, i.e. it belongs to an approximate (due to 
the non-uniform nature of the FMO complex) invariant sub- 
space whose concepts have been discussed here in the con- 
text of fully connected networks. Indeed, a numerical simu- 
lation shows that after a time t ~ 1 ps, the system rapidly 
evolves into a pure state where the main population is approx- 
imately evenly distributed in the sites 1 and 2 and where the 
amplitudes have a relative phase close to tt which in turn leads 
to destructive interference for transitions to level 3. It is this 
destructive interference that needs to be disrupted by a mod- 
erate level of dephasing noise on the sites 1 and 2. Indeed, 
for t = 5 ps, the noiseless evolution yields Psink = 0.566, 
while dephasing noise only on site 1 with 711 = 2.985 ps^^ 
yields Psink ~ 0.730, while dephasing noise only site 2 with 
722 = 20.111 ps-i yields = 0.772. 

Furthermore, it is crucial to note that the energy of site 3 
differs from that of sites 1 and 2 to an extent that the Hamilto- 
nian matrix elements between sites 1 and 2 are not sufficient 
to lead to strong transfer to site 3. This is why a moderate 
amount of dephasing on site 1, 2 or 3 may assist the transfer 
process using line broadening to strengthen this transfer chan- 
nel. Indeed, for < = 5 ps, dephasing only on site 3 alone, of 
733 = 35.623 ps~^, already yields Psink ~ 0.658. 

It should be noted that it is not desirable for the system 
to transfer population into the sites 4, 5, 6 and 7 as this will 
lengthen the time until eventual transfer into the reaction cen- 
ter via site 3. Hence strong dephasing on those sites would be 
ideal to suppress this undesired exploration of the sites 4, 5, 6 
and 7. As both sites 1 and 2 are most strongly coupled to level 
6, amongst these levels the strongest dephasing would be ex- 
pected for that site. For < = 5 ps, a dephasing on site 4 alone, 
with rate 744 = 20.551 ps^^, would yield Psink — 0.644, 
while dephasing on site 5 alone, with rate 755 = 32.076 ps^^, 
would yield Psmk = 0.718, and dephasing on site 6 alone, 
with rate 766 = 25.43 ps~^, would yield Psink = 0.626. 

We have also determined numerically the optimized de- 
phasing rates for t ~ 5 ps, where the optimal dephas- 
ing rates yield Psink ~ 0.903 as compared to Psink = 
0.566 in the dephasing free case. Those are given by 
{0.157, 9.432, 7.797, 9.432, 7.797, 0.922, 9.433} ps'K 

Finally, we would like to explore briefly to what extent non- 
local, i.e. correlated, noise, may lead to further enhancements 
of the transport performance of the system. To this end, we 
perform a numerical optimization for i = 5 ps, now mak- 
ing use of all the jmn in Eq. (14) rather than simply the di- 
agonal elements that describe local dephasing only. We find 
that with non-local dephasing noise the best transport proba- 
biUty into the reaction center is given by Psink = 0.931 as 
compared to Psink = 0.903 without correlations. This small 
additional transfer probability suggests that, while non-local 



effects may play some role for exciton transfer in the FMO 
complex, it is not expected to lead to decisive improvements 
in this setting. As a final remark, the basic mechanisms that 
we have elucidated in previous sections all suggest a weak or 
moderate dependency on the noise level where only vanishing 
or extremely strong dephasing noise lead to a deterioration 
of performance. For a broad intermediate level of noise and 
static disorder, the effects of suppression of destructive inter- 
ference and the enhanced coupling due to line broadening are 
active. This is corroborated here where we show that these 
optimal dephasing values are highly resilient under variations 
in the system parameters, like site energies, coupling strengths 
and dephasing rates. This is proved by the fact that enhanced 
transport in the FMO complex remains essentially unaffected 
in the presence of random variations (20%) in site and bond 
energies, as shown in Fig. (6). As already demonstrated in 
Fig. (3) for the fully connected network, a numerical study of 
the FMO complex also shows that variations in the dephas- 
ing rates around the optimal values do not affect the transfer 
probabilities significantly. This strongly suggests that the ex- 
perimental results recorded for samples at 77 K would also be 
observable at higher temperatures. 



B. Entanglement dynamics 

Let us consider the presence of entanglement in the FMO 
complex in the presence of dephasing noise. We study the 
dynamics of a single excitation, initially injected in site 1, 



1-(2-7) 

{1-2)-(3-7) 

{1-3)-(4-7) 

(1-4)-(5-7) 

(1-5)-(6-7) 

(1-6)-7 




0.4 0.6 
time (ps) 



FIG. 7: Entanglement, as quantified by the logarithmic negativ- 
ity [38], across different bipartitions, as specified in the inset, for 
optimized local dephasing rates at each site. The curves are for 
the 6 spUts of the form (1, • • ■ , fc)-vs-(fc + 1, • ■ ■ ,7), with k £ 
{1, ■ ■ ■ , 6}, of the FMO complex sites. Under the action of local 
dephasing, entanglement in our model persists up to 1 ps and there- 
fore quantum correlations are dissipated well before the excitation 
transfer (yellow line) is completed. 
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employing the master equation and Hamiltonian introduced 
in the previous subsections, in both the cases of local and 
spatially correlated dephasing. For the case of local dephas- 
ing, we use the optimized dephasing rates as above, while for 
spatially correlated case the diagonal elements of the matrix 
7m„ are again the optimized local dephasing rates and the off- 
diagonal elements are random numbers. Results for the local 
dephasing case are presented in Fig. (7). The initial exciton 
delocalization rapidly creates entanglement, quantified by the 
logarithmic negativity of the state, between site 1 and the rest 
of the complex (black line), which degrades as time elapses. 
Let us remind that the logarithmic negativity of a state p of 
a bipartite system AB is defined as log | \p^^ | |i, where Ta is 
the partial transpose operation with respect to the subsystem 
A and || • ||i denotes the trace norm There is also entan- 
glement across other bipartitions, with a profile that is consis- 
tent with the expectation of limited spatial distribution of the 
quantum correlations. The time evolution under spatially cor- 
related noise is broadly comparable to the uncorrected case 
shown in Fig. (7), although there is slightly more and longer 
lived entanglement in the correlated case. These results are 
in qualitative agreement with recent studies on entanglement 
dynamics in bipartite qubit systems subject to forms of corre- 
lated (non-Markovian) noise [39]. 

Notice, however, that the constraint on the matrix ^„m 
(which has to be positive semidefinite), i.e. 7 > 0, is essential 
here, while violation of this constraint leads to a time evolu- 
tion in which entanglement persists for much longer time but 
it does not preserve positivity and is hence not physical. This 
highlights the importance of considering only time evolutions 
that are completely positive and trace preserving and that even 
minor violations may lead to significant deviations in the en- 
tanglement dynamics. 

The present analysis was carried out in the single exciton 
sector in which the presence of coherence (off-diagonal ele- 
ments) is equivalent to the presence of entanglement. A more 
detailed analysis, relaxing this condition and studying the dy- 
namics under realistic laser excitation as well as driving from 
the antennae complexes in the presence of spatial and tempo- 
ral noise correlations will be presented elsewhere [40]. 



A. Destructive interference suppression 

The crucial importance of destructive interference between 
pathways and its suppression by dephasing noise can be illus- 
trated by a fundamental experiment (see upper part of Fig. (8)) 
that is capable of elucidating many of the basic aspects of 
quantum mechanics: the balanced Mach-Zehnder interferom- 
eter [41]. With a photon incident on the upper input port of 
a 50/50 beamsplitter, one finds that it will always emerge in 
the lower output port of the second beam-splitter while it is de- 
structive wave-interference that prevents the photon to emerge 
in the upper port. This destructive wave-interference may be 
perturbed in two fundamental ways, each of which serves to 
extract the which-path information about the photon [42, 43] 
- see lower part of Fig. (8). Firstly, local dephasing in one 
of the arms of the Mach-Zehnder interferometer (e.g. ran- 
dom fluctuations of the refractive index) will vary the phases 
and therefore inhibit the destructive interference. Secondly, 
when the photons in the beam-splitter suffer a path-dependent 
deterministic frequency, time or polarization shift, then wave- 
interference is inhibited because it is now possible to identify 
which path was taken. The essential insight gained by this 
simple example is the fact that decoherence and static disor- 
der, by inhibiting destructive interference, may have the effect 
of opening up additional pathways for propagation in the sys- 
tem (here the photon ending up either in the upper or the lower 
output port). 
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V. A TRANSPARENT PHYSICAL PICTURE 



Time 



The purpose of this section is to summarize the finding of 
previous sections and present, in a qualitative and transpar- 
ent manner, the mechanisms that may underlie the dephasing 
assisted transport of excitations in quantum mechanical mod- 
els of fully connected networks and, in particular, in light- 
harvesting complexes. 



FIG. 8: (a): A single photon entering a balanced Mach-Zehnder in- 
terferometer will always emerge in one output port as the path to the 
upper detector is blocked by destructive interference, (b): If the in- 
terferometer is unbalanced, as a result of a noisy process, the condi- 
tion for destructive interference is inhibited and photodetections are 
observed in both ports. Thus, noise may open additional paths for 
propagation. 
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Little overlap -> weak coupling 
Large overlap strong coupling 



Large gap weak coupling 
Small gap -> strong coupling 



FIG. 9: Left: Due to energy fluctuations, dephasing leads to a broad- 
ening of energy levels and hence increased overlap between sites. 
Right: Viewing these fluctuation dynamically, one finds that the en- 
ergy gap between levels varies in time. A non-linear dependence of 
the transfer rate on the energy gap may therefore lead to an enhance- 
ment of the average transfer rate in the presence of dephasing noise. 



broadened lines on the other. The fundamental character of 
these processes suggests that they play a role in a wide variety 
of transport processes whose physical and biological realisa- 
tions extend well beyond the specific example of photosyn- 
thetic complexes discussed here. Our results suggest that Na- 
ture may actively utilize the fundamental sensitivity of quan- 
tum systems to noise and that the way Nature seems to exploit 
environmental noise is intrinsically robust to changes in the 
properties of the dissipative quantum network. This robust- 
ness corroborates the idea of the general applicability of the 
described underlying principles in a variety of transport net- 
works and the possibility to exploit these noise-assisted pro- 
cesses for achieving robust and efficient energy transfer in ar- 
tificial nano-structures. 



B. Line broadening 

In addition to the inhibition of wave-like interference, de- 
phasing has another effect that may already be seen to be rele- 
vant on a classical level but which can also be expected to play 
a role in quantum systems. Indeed, dephasing may arise from 
the random fluctuation of the on-site energies which, in turn, 
implies that the energy gap between neighboring sites will 
fluctuate. As the effective coupling rate between sites usually 
depends non-linearly on the energy difference between neigh- 
boring levels, these fluctuations may indeed have the effect of 
enhancing the average transport rate between sites (rh.s. of 
Fig. (9)). Another way of expressing this observation is by 
saying that dephasing leads to line broadening so that excita- 
tion lines of interacting sites may overlap more strongly and 
hence enhance transport (l.h.s. of Fig. (9)). Hence, again, 
propagation paths that may not be open, or are weak, in the 
noise free case may be opened up (enhanced) in the presence 
of dephasing. 



VI. CONCLUDING REMARKS 

We have analysed transport assisted by dephasing, and 
more general types of noise, in dissipative networks such as 
those describing the transport of excitons in photosynthetic 
complexes [18, 19] and also considered the presence of en- 
tanglement. We have introduced an analytical technique that 
allows for the identification of invariant subspaces of network 
Hamiltonians and therefore the possibility of supporting noise 
assisted transport. We have further clarified the fundamen- 
tal mechanisms that lead to this phenomenon and provided 
simple examples to reveal their basic character. Of key im- 
portance are the suppression of destructive interference by 
dephasing noise, static disorder in energy levels and differ- 
ences in propagation speed along different paths on the one 
hand, and the enhanced overlap between sites with dephasing 
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APPENDIX A: COMPLETE SOLUTION OF THE MASTER 
EQUATION 

1. Calculating the population transfer 

We will now present a solution of the Markovian master 
equation for a fully connected network (FCN) with F, = 
r, 7i = 7, and uji = ui. In terms of the density matrix ele- 
ments in the site basis pij (t), the equations of motion are, 

Pij = - ["^T + rN+iiStN + Sjn) + 2j - 2^5ij] p,j 

N 

Poo = 2^Tjp.jj , 

where poo is the environmental population, the dots represent 
time derivatives, and the time labels, e.g. (i), have been sup- 
pressed in order to neaten the presentation of this and all sub- 
sequent equations of motion. We now introduce a set of col- 
lective variables given by 

N N N 

= Py and A = ^ i?, = ^ p,j , (A2) 

j = l i=l i,j=l 
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and reexpress the equations of motion in terms of these col- 
lective variables, i.e.. 

Pi, = -2Vp,, + U{R, -R,) , i^N 

= -2(r + 7)p,, + iJ{R, - Rj) , i^N.j^N 
pm = -{2T + 2-i + VN+i)PrN +^J{R^- Rn) , (A3) 

PNN = —2{TN+i+r)pNN + iJ{RN — R]y), 

Poo = 2r Trp . 

where Ri denotes the complex conjugate of Ri. These equa- 
tions are supplemented by the equations of motion for Ri and 
Rn, i.e., 

R, = -iJA + iJNR,-2{r + j)R,-TN+iP,N 

+2jpu, (A4) 

Rn = -iJA + iJNRN - {2r + 2-f + rN+i)RN 

+ {2-f^TN+i)pNN, (A5) 

and the equation of motion for A, i.e., 

A=-2(r + 7)A - Tn+1 {Rn + i?7v) + 27 Tip . (A6) 

The final term is the total population in the network, which is 
not preserved due to the losses to the environment and sink. 
However the following relation must always hold, i.e. 



1 = Trp + Poo + P. 



'sink 



(A7) 



Note now that the collective variables Ri will in general be 
complex, while A is always real. In particular, let 



Rn ^X + iY 



(A8) 



and this results in the following differential equations which 
form a closed system: 

A = -2{T + -f)A-2TN+lX + 2-f{l- poo-Ps^nk), 

X = -i2T + 2-/ + rN+i)X + i2j-TN+i)pNN-JNY, 
Y = -{2T + 2-f + rN+i)Y + JNX -JA, 
PNN = -2{T + Tn+i)pnn ~ 2JY, 

Poo = 2r(l - Poo - Psink), 
Psink — SFtv+iPattv ■ 



The initial conditions considered in the main text are 



A = l, X^O,Y = 0, pNN = 0, Poo = 0, p 



sink 







This system of coupled differential equations can be converted 
into a set of algebraic equations via the Laplace transform C. 
Decoupling these equations, one finds that the Laplace trans- 
forms of the dynamical variables contain only simple poles 
and the inverse Laplace transformation can be performed an- 
alytically. The analytical results presented in the paper are all 
derived from this procedure. In particular, the full dynamical 



problem can be reduced to the solution of the following set of 
equations for the Laplace s-domain variables A = £[A(i)] , 
i.e., 

(s + 2r + 27) A + 2rjv+iX + 27p,i„fc + 27P00 - 27/s -1 = 0, 
(s + 2r + 27 + rAr+i)X + (rjv+1 - 2'y)pMN + J NY = , 
(s -I- 2r + 27 + r,^+i)Y + JA- JNX = , 
(s + 2r + 2rjv+i)pjv]v + 2JY = , 

(s + 2r)poo + 2rp,i„fc - 2r/s = o , 

sPsink - 2r]v+i/5iv]v = . (A9) 
The quantity of eventual interest is Psink^ but the last equa- 
tion can be used to eliminate Psink in favour of pnn, leaving 
five linear equations involving an equal number of variables. 
When determining the poles of the Laplace transform of pat at, 
one finds that the characteristic equation is therefore one of 
fifth-order, and consequently its roots cannot be generally ex- 
pressed in analytical form. However, some special examples 
in which the dynamics of the full system can be analytically 
derived will be shown in the following. 

2. The full density matrix dynamics without losses or pure 
dephasing 

As an example of how the full dynamics and final state 
of the system can be extracted from the system of equations 
given above, we present an analysis of the case F = 0, 7 = 
0. The coupled equations for the collective variables in the 
Laplace domain given in Eqs. (A9) are reduced to 

sA + 2Tn+iX - 1 = 0, 

{s + Tn+i)X + TN+IPNN + JNY = 0, 

{s + Tn+i)Y + JA- JNX ^ Q , (AlO) 

{s + 2Tn+i)~Pnn + 2JY = 0, 

SPsink — 2Tn+iPNN ~ . 

Note that in the absence of dissipation to the sink, the system 
could be solved completely in terms of the collective variables 
via just the first three of the governing equations. The dissipa- 
tion to the sink is a local process and requires that we include 
Pn n in the system of equations in order to find the full dy- 
namics. Let us point out that introducing other local variables 
into the density matrix equations, e.g. pure dephasing, local 
dissipation, different on-site energies etc., would similarly re- 
quire extra local variables to enter into the collective equations 
of motion, e.g. the sink population and the environmental one. 
For the case considered here, one finds that the Laplace trans- 
form of PNN, is given by 

2j2(s + rjv+i) 



PNN 



{s - st){s - Si )(s - 4)(s - S2 ) ' 
where the are the four roots of the equation 

= s(.s + Fw+i)'(s + 2FAr+i) 

+ J^iNs + 2rN+i)iNs-2TN+i+2NrN+i) 



(All) 
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The laplace transform has four simple poles, and the inversion 
is now straightforward. The solution in the time domain is 
given in terms of the residues of as 



Pnn[ ^ 

i=1.2 

(A12) 

The roots above can be determined analytically, and can be 
written in the form. 



Si = -rjv+i ± A+ , 



(A13) 



where 



Note that pNNit) vanishes as t ^ oo, as one would expect. 
This is ensured by the relation F^v+i > A±, which using 
Eqs. (A 13) can be shown to be always valid. Calculating the 
residues of these poles, one finally obtains pNN{t) from Eq. 
(A12), 

2j2e-r"+i* [cosh(A+t) - cosh(A_t)] 



PNN(t) 

(A14) 

Once pNN{i) is determined, psink is obtained easily via. 



Psink 

(t) = 2rjv+i / pNiq{t')dt' . (A15) 
Jo 

Taking the limit t oo, one recovers the result found by both 
numerical simulations and the invariant subspace approach 
(see below), which is 



Psmfc(oo) = 



1 



N - 1 



Having determined pNN{t) we can now, in theory, deter- 
mine the rest of the density matrix elements. Formally one 
has to substitute the solution for pNN{t) back into the system 
of equations in order to determine the collective variables A{t) 
and i?7v (t)- Doing this explicitly is rather tedious, and we will 
focus instead on some general features of the solutions. First, 
consider the equation of motions for Ri{t) and pi^ {t), 

i?j = iJNRi-TN+iPiN - iJJ^ , 

PiN ~ ~^N+lPiN + iJ{Ri — Rn) ■ 

These can be combined to obtain a closed second-order dif- 
ferential equation for pay, 

PiN + {Tn+1 - iJN)piN - iJTpf+iiN - l)piN 

= J2(A - NRn) - URn ■ (A16) 



The initial conditions are 

p,Ar(0) = 0, P,n{0)=0, /9ijv(0)=0, piN{0)=iJ 

Looking at Eq. (A16), one can see that the rh.s. is a fixed 
function of the previously determined collective variables i?jv 
and A. Given that for i 7^ 1 all the initial conditions are iden- 
tical, all PiN must be identical functions of time. However, 
the different initial conditions for pn mean that this coher- 
ence can be different from the others. A similar second-order 
equation for Ri can also be obtained, showing that all Ri are 
the same except for Ri and Rpf. The equations of motion for 
the remaining density matrix elements are. 



Pij = iJ{Ri - Rj), 



(A17) 



which implies that all matrix elements are identical functions 
of time for i,j ^ N, except pn and pij. A further implica- 
tion is that pij ~ pik for all j,k ^ 1, N. With these results 
one can obtain several results concerning the stead-state pop- 
ulations and coherences. The steady states of the collective 
variable equations of motion lead to the final values, 

a:==o, y = o, PNN^o, a = o. (ais) 

Substituting these results into Eqs. (A16), one finds that all 
PiN = in the steady-state. The steady-state populations 
must also obey the conservation of population sum rule. 



N 



1 = ^ P]j + Psink , 



(A19) 



which using the results pnn = 0, Psink = {N — 1) ^, and 
Pa — Pu for i,l ^ 1,N implies the condition. 



iV-2 
N - 1 



= iN- 2)pu + Pii 



(A20) 



Setting the time derivative pn = 0, the following steady-state 
condition must hold. 



{N^l)pu = pii + (iV-l)p,i 



(A21) 



Setting Ri ~ one also finds that Ri = in the steady-state. 
This imphes. 



Pij 



Pii 
' N -2' 



{All) 



Combining these relationships, the final populations and co- 
herences are 



Pii 



A^-2 
N - 1 



Pii = 



N -2 

(iV- 1)2' 



Pij — Pii ' 



[N-iy 
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3. The full density matrix dynamics with losses and pure 
dephasing 

In the most general case the equations of motion for the in- 
dividual matrix elements of the density matrix can be written 
as 



Pi] 



Poo 



- [Ti + Tj + TN+i{StN + Sjn) + li + 7i - 27i(5y 
- cjj)] /9y Pa - X! P^i ) ' (^23) 



N 



(A24) 



Under the assumption that all the sites have equal energies, de- 
phasing and dissipation, Eqs. (A9) can be solved completely 
in terms of a few collective variables but analytically only in 
the Laplace domain. Indeed, the sink population is then ex- 
pressed, in terms of the Laplace variable s, as 



Psink{s) = ij'^TN+l 



sA(s) 



(A25) 



with 



A(s) = (s + 2r)(s + r,)(s + rb)2(.s + r,) 

+ 4j2rAr+i(47(s + 7 + 2r) - J^TN+iis - 27 + 2r) 



J^N\s + 2T){s + Ta){s + r,) , 



(A26) 



where 



r„ = 27+2r, Lfc = 27+2r+rAr+i, r, = 2r+2r;v+i . 

(A27) 

We can now use the above equations to solve for the sink pop- 
ulation at any time, but that involves solving for the roots of 
A(s) — 0, which being a fifth order polynomial cannot be 
solved in closed form. Here, we stress on the asymptotic value 
of the sink population which is easily obtained via the final 
value theorem and is given by 



r„r. 



Psink{t = OO) = \iin[s Ps.nk] = 4:J^'^N+1^, 



A(0) 



(A28) 



This expression is considerably general, and in particular one 
can recovers the noiseless case, i.e. 



Psink{oo)= lim 

7^0, L 

r 



4J i N+1 



A(0) 



1 



N -1 



(A29) 



For the case of purely local dephasing, i.e. F = 0, one finds 

= 1 . (A30) 



Psmfc(oo) = lim AJ'^Tn^"'^'' 



r^o 



A(0) 



4. An invariant subspace approach for FCNs 

We have shown that the full density matrix equations can 
be solved exactly in a number of cases, and have used these 
to obtain - somewhat laboriously - several of the results found 
numerically. This section presents a more detailed develop- 
ment of the invariant subspace method, which, as discussed in 
the main text, is a powerful method that allows one to obtain 
the results given above in an intuitive way and which can be 
generalised to deal with situations that cannot be easily tack- 
led by solving the dynamical equations. 

For the case of Tj = 0, 7j = 0, the initial state of the sys- 
tem can be expanded in terms of invariant and non-invariant 
defined in the main text. This expansion is explicitly given by 



|1) 



N - 1 



Invariant 



1+) 



^(1^,) I +VN\^)-\N) 

j=2 



(A31) 

The invariant states are eigenstates of the Hamiltonian and are 
decoupled from the dissipative dynamics that acts at site N. 
The only non-trivial time development comes from the term 
in the expansion that lies outside of the invariant subspace. 
Defining the state | +) as. 



|+) = ViV|0)-|7V), 

one can show that the amplitudes of a state of the form 
\ip){t) = a{t)\+) + b{t)\N) obey the following closed sys- 
tem of equations of motion. 



d(t) 

m 



JiN- 
J{N- 



2) J 
1) 



ait) 

m 



(A32) 



These equations show that spectral weight oscillates between 
1+) and \N), and will thus be entirely lost to the sink as t — > 
oo due to the irreversible relaxation which acts on site A^. As 
all the weight held in the non-invariant subspace is initially 
held in |+), we can immediately predict from the initial state 
expansion that. 



Ps 



c(oo) 



N - 1 



1 



(TV -1)2 (iV-l)2 iV-1 



(A33) 



The arguments given above can also be adapted to analyze 
the effects of local perturbations on the sites. These pertur- 
bations may be changes to the site energies or changes in the 
hopping rate between the sites and a given site. The change to 
the hopping rate of site (3 can be described by the new Hamil- 
tonian, 



J i' 



(A34) 
(A35) 
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We now introduce the reduced invariant subspace of the per- 
turbed Hamiltonian as discussed in the main text. For a net- 
work with D perturbed sites described by states \d), the initial 
state can be expanded in the invariant states as. 



ID 



1 



N ~D-l 



Y,\^,) + ^\<j^)-\N)-Y^\d) 



D 



(A36) 

where the tilde on the sum denotes that all which con- 
tain weight on a perturbed site are excluded. Once again one 
can show that the Hamiltonian dynamics leads to oscillations 
of the spectral weight initially held in non-invariant part of 
the expansion of |1) amongst the states \d), andjiV). 
The leakage to the sink through site therefore causes the 
eventual transfer of all of this weight into the sink. The fi- 
nal weight transferred to the sink is therefore Psink{oo) = 
{N - D - 1)-'^{D\D), which gives 

Psinfc(oo) = ^ _ ^ _ _^ . (A37) 

The approach outlined above, which essentially consists of 
writing the initial condition in terms of degenerate, invari- 
ant (to dissipation) states and a remaining non-stationary part, 
fails if the perturbed site is the injection site. The reason is 
that in this case one can no longer write the initial condition 
as we did before. If one tries to expand the initial state in the 
invariant states, one finds that 



ID 



N\4>) 



(A38) 



Comparing with Eq. (A36), one can see that there is no invari- 
ant subspace, the wave function only contains a non-stationary 
component which has precisely the form we have argued al- 
ways decays to zero. Therefore one finds that all the excitation 
is transferred to the sink and psink = 1 for any local perturba- 
tion on the injection site 1. 

5. Invariant subspaces for general networks 

Consider a network of N sites with a general Hamilto- 
nian that has D degenerate orthonormal eigenstates 1^"^) and 
A^—Z? non-degenerate orthonormaleigenstates |$fc). We shall 
now present a method to construct a generalised invariant sub- 
space that is unaffected by dissipation at site N (we assume 
no other sources of noise or losses in the system). Within the 
degenerate and non-degenerate manifold of states, isolate all 
states with no overlap with \N). These states are by definition 
in the invariant subspace. We now take linear combinations 
of the remaining degenerate eigenstates to create a new basis 
using the transformation 



l*() = 



D' 

E 



Ulk\^k) , 



where k runs over all D' degenerate eigenvectors with finite 
(AjV'fc). The overlap of the new basis with \N) is given by. 



1 ^' 



(A40) 



fe=i 



There are D' states I^P/) where the index / = 1 . . . D'. For 



I < D' the overlap is zero as ^ 



D' 

fe=i' 



0, and therefore 



this transformation creates D' — \ states in the generalised 
invariant subspace with (Aji/v) = 0. The transformation 
also generates one state that is outside the subspace given by 

\^D') = -D'"^ Efli(^l*fe>"^l*fc>- Thus an invariant sub- 
space can be generated for any network with degeneracy. This 
procedure systematically generates the invariant subspace, al- 
though it should be noted that basis is not unique as various 
other linear combinations of the degenerate eigenstates can 
be made. This is clear if one follows this procedure for the 
FCN, as the states generated here are not the pair eigenstates 
but rather delocalised states. We have also ignored the case of 
multiple degenerate manifolds, although the generalisation of 
the procedure to this case simply involves performing trans- 
formations of the form given in Eq. (A39) to each manifold. 

In order to study the transport dynamics of the network, we 
express the initial condition of the system 1 1) by an expansion 
in the original orthonormal eigenstates of the system. We then 
transform the degenerate eigenstates to the new eigenbasis us- 
ing the inverse transformation of Eq. (A39) and group the 
states into the generalised invariant and non-invariant groups. 
The part of the initial condition that can be represented by 
elements of the invariant subspace is unaffected by the dis- 
sipation at state |A) and the initial spectral weight will be 
preserved in the network. Using the transformation above, the 
amount of the initial excitation that is transferred to the bath is 
given by the weight held outside of the generalised invariant 
subspace. The population transferred to the sink is therefore. 



D \ D 

-2 



\k=l 



k=l 



Ps^nk = [ J2{N\^,)\-^ I \J2{N\^k){^k\l) 

N-D 



fc=i 



where it is assumed that there are no naturally occurring in- 
variant eigenstate, i.e. D' = D. 
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